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Excitation spectra of S — 1/2 and S = 1 frustrated Heisenberg antiferromagnetic chains with 
bond alternation (explicit dimerization) are studied using a combination of analytical and numerical 
methods. The system undergoes a dimerization transition at a critical bond alternation parameter 
S = S c , where <5 C = for the S =1/2 chain. The SU(2)-symmetric sine-Gordon theory is known to be 
an effective field theory of the system except at the transition point. The sine-Gordon theory has a 
SU(2)-triplet and a SU(2)-singlet of elementary excitation, and the mass ratio r of the singlet to the 
triplet is \/3- However, our numerical calculation with the infinite time-evolving block decimation 
method shows that r depends on the frustration (next-nearest-neighbor coupling) and is generally 
different from \f%. This can be understood as an effect of marginal perturbation to the sine-Gordon 
theory. In fact, at the critical frustration separating the second-order and first-order dimerization 
transitions, the marginal operator vanishes and r = \/3 holds. We derive the mass ratio r analytically 
using form-factor perturbation theory combined with a renormalization-group analysis. Our formula 
agrees well with the numerical results, confirming the theoretical picture. The present theory also 
implies that, even in the presence of a marginally irrelevant operator, the mass ratio approaches V3 
in the very vicinity of the second-order dimerization critical point S ~ 5 C . However, such a region is 
extremely small and would be difficult to observe numerically. 
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I. INTRODUCTION 

Techniques of field theory have achieved growing suc- 
cess in interpreting physical properties in low dimensional 
magnets. The achievement stems from the close connec- 
tion between one-dimensional quantum spin models and 
their effective theories. In particular, 5=1/2 Heisen- 
berg antiferromagnetic (HAF) chains with various per- 
turbations are important and also relevant for experimen- 
tal studies of one-dimensional magnets. The bosoniza- 
tion schemeEl is useful for analyzing these systems. A 
HAF chain with bond alternation, or under a staggered 
field, is described effectively by the sine-Gordon (SG) 
field theory. Elementary excitations in these systems are 
a soliton, an antisoliton, and breathers (bound states of 
the seiijton and antisoliton). Materials such as Cu ben- 
zoateBB and KCuGaFgP are described by HAF in a stag- 
gered field, and the soliton gap calculated from SG the- 
ory explains well the experimental results. For dimerized 
chains, the gap formula as a function .of dimerization S 
with logarithmic correction is obtainecB: <5 2 / 3 /| logi^ 1 / 2 , 
or it can also be represented as an effective power-law 
form, with a renormalized exponent which deviates from 
2/3. a Refined logarithmic correction is given in Ref. @. 
Dimerized spin chains are an appropriate model for spin- 
Peierls materials such as CuGe03@ or Ni compounds. El 

There are also a number of numerical studies on the 
frustrated HAF chain with next-nearest-neighbor cou- 
pling. We consider the Hamiltonian 

n = J E K 1 + (-tftysi • + aS i ■ s i+*\ > (!) 
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where J > 0. The next-nearest-neighbor coupling a > 



introduces frustration. 

This model exhibits a dimerization transition at 8 = 5 C . 
For 5 = 1/2, the transition point is always 5 C = 0, since 
the Lieb-Schultz-Mattis theorem implies either gapless 
excitations or two-fold degeneracy of the ground states 
at S = 0. In fact, on the undimerized line 6 — 0, there, 
exists a critical frustration parameter a c ~ 0.2411.E2IEJ 
For a < a c the system is a gapless Tomonaga-Luttinger 
Liquid (TLL); that is, the dimerization transition at 
5 = S c = is of second order. In contrast, for a > a c , the 
ground state is doubly degenerate, exhibiting a sponta- 
neous dimerization. This implies a first-order dimeriza- 
tion transition at S = S c = 0. 

For 5 = 1, on the other hand, 5 — (for a small 
a) belongs to the Haldane phase and does not repre- 
sent a transition line. Instead, a dimerization transition 
betwepp-the Haldane phase and the dimerized phase oc- 
cursEjTJ at a finite S c , which depends on the frustration 
a. Although the shape of the phase diagram is thus dif- 
ferent, the topology of the phase diagram is rather similar 
to that for S — 1/2. In fact, also for 5 = 1, there is a 
critical frustration a c ; the transition is second order with 
the critical point described by a TLL for a < a c , and first 
order for a > a c . 

In the neighborhood of the gapless TLL line, the 
system acquires a small excitation gap, and would be 
described by the SG theory. Since our model (|l|) is 
SU(2)-invariant, the SG theory should also have SU(2)- 
symmetry. As a consequence, the mass ratio r of the 
second lowest breather to the soliton should be \Z3. 

However, numerical results for 5 = 1/2 chainsll3 show 
that r generally does not agree with the SG theory pre- 
diction y/3. While r depends only weakly on 5, it does 
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vary as a function of a. Only near the critical frustra- 
tion a — a c does r agree with the SG prediction v3- 
In Ref. [l^, it was pointed out that a marginal operator 
exists as a perturbation to the SG theory, and it would 
shift r from y3- However, how exactly the mass ratio r 
is affected by the marginal operator was not clarified. 

The effect of the marginal perturbation to the SG the- 
ory on the mass ratio was discussed in terms of form- 
factor perturbation theory (FFPT) in Ref. [I?]. However, 
the theoretical prediction has not been tested. The mass 
ratio in the S — 1 case has also never been studied nu- 
merically. 

In this paper, we study numerically the mass ratio of 
elementary excitations and the ground phase diagram of 
the frustrated HAF with bond alternation (Q) for both 
S = 1/2 and S = 1. We employ the recently devel- 
oped infinite time-evolving block decimation (iTEBD) 
method, £3 which allows high-precision calculation of in- 
finitely long chains. The masses of elementary excita- 
tions are obtained from the asymptotic behavior of equal- 
time correlation functions, instead of extrapolation of the 
finite-size energy spectrum. We confirm previous results 
when they are available, and we obtain the mass ratio 
for S — 1 as a new result. Furthermore, we derive an 
explicit formula for the mass ratio r as a function of 5 
and a, by combining FFPT and renormalization-group 
analysis. This agrees well with the numerical results for 
both S — 1/2 and S = 1. Thus both cases are under- 
stood in terms of the unified framework of the SG theory 
with a marginal perturbation. 

This paper is organized as follows. In Sees, [n] and [II 
respectively, we review direct bosonization of the S = 1/2 
chain and derivation of the SG theory for general S case 
via the 0(3) nonlinear sigma model (NLSM). In Sees. |v| 
and 0, we present numerical study on the mass ratio and 
phase diagram, respectively for S = 1/2 and 5 = 1. We 
then discuss the mass ratio analytically based on FFPT 
and compare the theor etica l formula with the numerical 
results in Sec. VI. Sec. VII is devoted to conclusions. 
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(2na) 



dx cos (20) 



(2ira) 



dx cos(40). 



(2) 



Irrelevant terms are omitted here, u and K denote spinon 
velocity and the Luttinger parameter, respectively. At 
the SU(2)-symmetric Heisenberg point, u = ira/2 and 
K = 1/2. Since the operator e 19 *^ has scaling dimen- 
sion Kq 2 /4, the cos(20) term is relevant while the cos(40) 
term becomes marginal. The 51 term arises from the 
bond alternation (i.e. dimerization) . 32 is known to de- 
crease with increasing a and vanish at a — a c where the 
transition from TLL to the self-dimerized phase happens. 
Thus, coupling constants g\ and g^ are proportional to 5 
and a — a c , respectively. When 51 ^ and 52 = 0, (||) is 
equivalent to the SG model. It is an exactUwplved model, 
and the excitation spectrum is obtained£3'E2l There ap- 
pear three types of elementary particles: a soliton, a cor- 
responding antisoliton, and breathers. The number of 
breathers is [2/ K — I], where [A] stands for the integer 
part of A. The mass of soliton Ms and the 71-th lightest 
breather Mb„ are related through the formula 



M B „ = 2M S sin 



A/K -2 



n = 1, 



(3) 

According to (j^), in HAF chain with dimerization (K — 
1/2), the soliton, the antisoliton and the first breather 
form triplet while the second breather is a singlet which 
has v^-times as large mass as the triplet. Although 
the degeneracy of the triplet is protected thanks to 
SU(2)-symmetry, the mass ratio of singlet to triplet r = 
Mb 2 /Ms is subject to correction caused by the marginal 
term gi. 



III. SG THEORY VIA NONLINEAR SIGMA 
MODEL 



II. BOSONIZATION 

We first review the bosonization of a spin-I/2 chain. 
Spin operators are represented as 

S z . = -d x <t> + ax(-l) j cos(20) + • • • 
S+ =c w [b (-iy+b 1 cos(20) + -..] , 

where dual boson fields <j>, 6 satisfy the commutation re- 
lation [4>(x), 9(x')] = —iiri9(x — x') {d{x — x') is the step 
function) with x = ja (a is lattice spacing), ip and 9 have 
periodicity </j~</j + 7t, + 2ir. Effective Hamiltonian 
of XXZ chain with dimerization is written with cf> and 9 



S > 1/2 chains may be bosonized by introducing 
Hund coupling to 2S chains of spin-I/2. Each chain is 
bosonized separately, resulting in a theory of interacting 
2S boson fields£H In the low-energy limit, however, one 
of the linear combinations of the boson fields becomes 
important. The SG theory (or TLL) would emerge as an 
effective theory of this linear combination. 

However, it is rather cumbersome to pursue this ex- 
plicitly. As an alternative, the SG theory can also be 
derived from the 0(3) non- linear sigma model (NLSM). 
The 0(3) NLSM was derived in the semi-classical, large- 
S limit of the HAF chain. Nevertheless, it proved to be 
a useful effective theory even for S = 1. 

Let us define fields n{x) and l(x) by Sj/S ~ 
{-l)' J n(x) + l(x). Then the spin-5 HAF chain with 
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bond alternation (|T|) can be generally mapped to the 0(3) 
NLSM 

A ° = Ygl dTdT [ V ( d * n ^ + + 

where g = 2/5 is some coupling constant and w = 2J5 
is the spin-wave velocity. T = J drdxn ■ d x n x d T n 
represents the integer-valued topological charge and 9 = 
2itS(1 + 5). For the moment, let us assume that there is 
no frustration, a — 0. 

0(3) NLSM is known to be integrableB@ at 9 = 
and 7r. At 6 = mod 27r, the excitation consists of 
a triplet of massive particles. In contrast, the theory is 
massless at 9 = 7r mod 27r and the infrared fixed point is 
a SU(2)i Wess-Zumino-Witten model, a conformal field 
theory (CFT) with central charge c = 1. This is merely 
the TLL at the SU(2)-symmetric point K = 1/2. 

When bond alternation is absent (8 = 0), the system 
is massless (9 = it) if 5 is a half-odd-integer, while it is 
massive (0 = 0) if 5js an integer. This is the celebrated 
Haldane conjecture ,E-3 which is now established by inten- 
sive analytical, numerical, and experimental studies. 

It is also interesting to consider the effect of bond al- 
ternation 6. By changing 8 from —1 to 1, namely from 
the completely dimerized limit to the opposite completely 
dimerized limit, 9 passes the critical point, n mod 2-7T, 
2S times. Thus, on — 1_ .< 6 < 1, there are 2S suc- 
cessive phase transitions .EJ This could be understood as 
successive spontaneous breaking and restoration of hid- 
den symmetryo or more generally, symmetry-protected 
topological phase transitionsrjO 

For 5 = 1/2, the transition occurs only at 6 = 0, 
consistently with the direct bosonization analysis. For 
5 = 1, there are two transitions which separate the Hal- 
dane phase around 6 = from the dimerized phases. 
The critical points are, according to the above argument, 
given by S = ±S C = ±1/2. However, in reality, the loca- 
tion of the critical points is renormalized. It was shownta 
numerically that S c ~ 0.25 J. 

As discussed above, the critical point is described by 
the SU(2)-symmetric TLL with K = 1/2. By considering 
the possible perturbations to the TLL, the effective the- 
ory near the critical point <5 = <5 C is determinedEZI to be 
the SG theory with marginal perturbation (||) , which was 
derived previously for S = 1/2 by direct bosonization. 
Thus, the same theory (||) should describe the neighbor- 
hood of dimerization transitions for any 5. In the follow- 
ing, we shall investigate the systems with 5 = 1/2 and 
5=1 numerically, and verify this universality. 



IV. MASS RATIO AND PHASE DIAGRAM FOR 

5=1/2 



extract the excitation spectrum from the equal-time cor- 
relation function obtained by iTEBD, shown as follows. 

A single-particle excitation in the SG model can be 
parameterized by the rapidity 9, which defines its energy 
and wave number as Mq cosh# and (Mq/u) sinh#, respec- 
tively (Mo is the mass of the particle). The one-particle 
form factor of operator O is specified by 9 and the kind 
of particle a as Fo(9,a) = (0\O\9,a). O represents an 
operator which creates the single soliton, the antisoliton, 
or the breather. We can calculate the equal-time corre- 
lation function (O(r)O(O)) - (O(r))(O(0)) by inserting 
the resolution of the identity 1 = X^^Lo ^ n wnere Pn is 
the projection operator defined as Pq = |0)(0| and P n = 

h J2au-,a n I ^J)"' |0i» a i; • • • ;9 n ,a n )(9i,a 1 ; - ■ ■ ;0„,a„| 
(n > ljh-. Then, the leading order of the correlation 
functioned is 

<O(r)O(0)> - <O(r))<O(0)> 

w / ^e iM °'' sinhe /"|F o (0,a)| 2 . 

J 47T 

In the limit of / — > oo, it is calculated to beEl 



(O(l)O(0)) - (O(l))(O(0)) = (A(-l)' +B)- 



-i/t 



(4) 



consisting of a staggered and uniform part. We suppose 
that the effect of the marginal cos(4</>) term is renormal- 
ized into mass Mq and constants A, B. 

In this way, the mass can be extracted from the cor- 
relation function, which we calculate with the iTEBD 
method. The truncation dimension, the number of con- 
served states in evolution, is fixed to be 200, large enough 
for the iTEBD calculation in gapped systems. (S^Sf), 
(£%£%), (S$Sf) and ((S • Si)(S t ■ S l+1 )) are fitted with 
Ce~ 1 ^ /Vl for sufficiently large and even I. C(= A + B) 
and £ are fitting parameters. Then we can obtain the 
mass of the soliton, the antisoliton, and the first and sec- 
ond breathers, respectively, through the relation M = 
u/£. Note that M is a renormalized mass. While the 
value u for a = is obtained exactly from the Bethe 
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We study the excitation spectrum of the system nu- 
merically, and we focus in particular on the change of r 
due to the marginal term. We adopt a new strategy to 



FIG. 1. (Color online) Correlation functions (Sf)Sf) and 
(SgSf) calculated with the iTEBD method. The solid line 
represents the fitting with the function Ce~ 1 ^ / 
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FIG. 2. (Color online) (a) Triplet-singlet mass ratio r as a 
function of a and 5 in the S = 1/2 bond-alternating chain 
with frustration, (b) Phase diagram of the S — 1/2 bond- 
alternating chain with frustration. Solid and dashed lines 
represent second order (TLL, c = 1 CFT) and first order tran- 
sition, respectively. Universality class of transition changes at 
a ~ 0.25, where r becomes v3- 



ansatz, it cannot be for a^0. Yet, the value of u is not 
needed to calculate a mass ratio. 



Since St 



Y]j Sj commutes with the Hamiltonian, 



tot is a good quantum number. The ground state is in 
Stot = Hilbert space, and the soliton (antisoliton) is 
an excitation to the lowest energy level in 5 t z ot = 1(— 1) 
Hilbert space. Hence, their mass corresponds to the in- 
verse correlation length of the operator changing S* ot 
by ±lE2, i.e. (SfiSf) = (S%Sf). On the other hand, 
the first breather is the lowest excitation in S* ot = 
Hilbert space; it corresponds to (S^Sf). In the case of 
the antiferromagnetic XXZ model, the mass of the soli- 
ton / antisoliton and the first breather is different. For 
a Heisenberg chain, however, SU(2)-symmetry requires 
{S%S?) = (S%S?) = (5 2 5f), which indicates that the 
mass of the soliton, the antisoliton, and the first breather 
is all the same, and these three particles constitute a 
triplet. The second breather has to be a singlet, and the 
operator corresponding to it does not change S* ot . The 
most relevant operator with such properties is Sj ■ Sj + \, 
and we expect that the second breather corresponds to 
((S -S 1 )(S l -S l+1 )). 

An example of fitting for correlation functions is shown 
in Fig. @, (SgS?) and (SgS? ) calculated with the iTEBD 
method are equal up to eight digits, which is consis- 
tent with the SU(2)-symmetry. The solid line represents 
the fitting with the function Ce~ 1 ^ /VI. The correlation 
functions are well fitted with the function. 

We show numerically the calculated mass ratio r as a 
function of a and 5 in Fig. ||(a) . r is larger than 2 for a — 
(nonfrustrated HAF chain with bond alternation) and 
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FIG. 3. (Color online) (a) Triplet-singlet mass ratio r as a 
function of a and S in the S = 1 bond-alternating chain with 
frustration. The transition point S c corresponds to circles in 
(b). (b) Phase diagram of the S = 1 bond-alternating chain 
with frustration. The solid and dashed lines, which represent 
the second order (TLL, c = 1 CFT) and first order transition, 
respectively, are guides for the eye. The circles show tran- 
sition points determined from Fig. ^J. The universality class 
of the transition changes approximately at a = 0.3, where r 
becomes v3. 



decreases with increasing a. r becomes V3 at a ~ 0.25. 
It is very close to a = 0.2411, where the transition from 
TLL to the self-dimerized phase happens without bond, 
alternation, and the marginal cos(40) term vanishes.Eil 
This result indicates that the deviation of r from y3 is 
attributed to the effect of the marginal term. While r is 
subject to correction as a moves away from this point, 
its 5 dependence is quite small. 

A similar result was obtained through a gap evaluation 
by exact diagonalization.t£l However, the mechanism of 
the variation of r has not been made clear. We will the- 
oretically analyze the dependence of r on the frustration 
a later in Sec. 



VI 



The a- 6 phase diagram is shown in 
Fig. |2](b). Note that the universality class of the tran- 
sition from positive to negative 6 is oh-e, = 1 CFT for 
a < 0.25 and of first order for a > 0.25.1 



V. MASS RATIO AND PHASE DIAGRAM FOR 

S= 1 

Next, we numerically investigate the excitation spec- 
trum and the phase diagram of the S = 1 HAF chain 
with dimerization and frustration. The method for eval- 
uating particle mass is the same as for the 5 = 1/2 chain. 
As can be seen in Fig. [| when a is small enough, r is 
always larger than 2 at least in \S — S c \ > 0.005 (5 C can 
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FIG. 4. (Color online) (a) Correlation length £ as a function 
of 8 — 8 C and a. £ diverges when the transition is second 
order, (b) Dimerization order parameter \{(—iySj ■ Sj+i}\ as 
a function of 8 — S c and a. 
the transition is first order. 

from the divergence of £ or the jump of |((— l) 3 Sj ■ Sj+i)\. 



\{(—iySj ■ Sj+i}\ jumps when 
Therefore, 8 C can be determined 



be determined from the divergence of £ or the jump of 
|((— lp'S'j • Sj + i) \ as explained in the following part. See 
Fig. |j), and does not depend much on 8. Since particles 
heavier than 2Ms become resonance, the second breather 
cannot be a stable particle even in the vicinity of <5 C . The 
above result again seems inconsistent with the prediction 
in Ref . [HJ 

The deviation of r from V?> would be attributed to the 
existence of the marginal term as in the spin- 1/2 chain. 
We introduce the next-nearest-neighbor coupling a in or- 
der to confirm it. As shown in Fig. ^(a), r decreases 
with increasing a and becomes Vo around a — 0.3. The 
transition point 5 C from the Haldane phase to the dimer- 
ized phase also decreases, which is natural because next- 
nearest-neighbor coupling favors the dimerized phase. 
Fig. [| shows the behavior of correlation length £ and 
dimerization order parameter \((—iySj ■ Sj+\}\ near S c . 
£ diverges at S c for a < 0.3, which is not the case for 
a > 0.3. In addition, |((— iySj ■ Sj+±)\ jumps at 5 C for 
a > 0.3 while the variation is continuous for a < 0.3. 
These results indicate that the universality class of tran- 
sition at 8 C changes from c = 1 CFT to first order when a 
goes beyond 0.3. From the viewpoint of field theory, the 
cos (4</>) term changes from a marginally irrelevant to a 
marginally relevant operator at this point. The situation 
is very analogous to the spin- 1/2 case. The a- 5 phase 
diagram is summarized in Fig. 0(b). It is consistent with 
Ref. |l3|. The transitions along the lines 8 = and a = 
are studied in Refs. |l2| and [Hi], respectively. 



VI. MASS RATIO FROM THE FORM-FACTOR 
PERTURBATION THEORY 

Now let us discuss the variation of the mass ratio r 
theoretically. In Ref. [l7| the 8 dependence of r was dis- 
cussed as follows. The excitation structure at the very 
vicinity of S = 8 C would be described by the pure SG the- 
ory without the marginal perturbation; r is then equal to 
V3. On the other hand, 0(3) NLSM with 6 = also has 
triplet lowest excitation, which is smoothly connected to 
the triplet in the SG model thanks to SU(2)-symmetry, 
but does not have the second breather. Therefore, r in- 
creases as 8 decreases from 8 C to 0, and it exceeds 2 at 
some point. This argument was further augmented by a 
FFPT calculation in terms of-the marginal perturbation. 

However, their predictionsO do not seem to be consis- 
tent with numerical results. In the absence of frustration 
a, r is substantially larger than y/3 even when 8 is clos- 
est to 8 C within the precision of the numerical calcula- 
tions. This already contradicts the picture presented in 
Ref. 17, Moreover, the effect of the frustration a was not 



discussed. 

Here, we will improve the FFPT by supplementing it 
with a RG analysis. Let us define a dimensionless cou- 
pling constant y 2 = gi/^u). With the FFPT of the 
marginal operator in the SG theory, mass corrections 
arising from the marginal term y 2 to the triplet and the 
singlet, which we denote, respectively, as AM t and AM,, 
were foundO to be 



AM t = ±V3y 2 , 
AM, = l2V3y 2 . 



(5) 



Here, we argue that the renormalized coupling constant 
should be used for y 2 . In the following, we derive the 
renormalized form of y 2 . Since the system has SU(2)- 
symmetry, y 2 is renormalized according to the Kostcrlitz- 
Thouless renormalization equationtrEj, 



dy 2 
ds 



= yl 



(6) 



The solution of @ is y 2 = — l/(s + Const.). y 2 becomes 
a function of energy scale by the parametrization s = 
hx(E/A) (A is the infrared cutoff) as follows 

V2(E) = 



]n(A'/E)' 

Constant A' can be fixed from the condition that bare y 2 
corresponds to the original spin chain, where the energy 
scale is of order of J, i.e., y 2 {E ~ J) = Ci(a c — a), where 
C\ is a non-universal positive constant. Therefore, the 
renormalized form of y 2 is 



V2(E) = 



1 



ln(J/E) 



d(a c -a) 



When the system is renormalized until the energy scale 
is equal to the soliton mass, y 2 becomes y 2 {M$). From 
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irrelevant operator is renormalized to zero very slowly 
(logarithmically), and thus the mass scale must be expo- 
nentially small in order to probe this regime. This can 
indeed be seen in the logarithmic dependence of r on the 
soliton mass Mg in eq. (|t]). Thus, for a < a c , the mass 
ratio deviates very quickly from r — as S is shifted 
from the critical point 5 C . As a consequence, it would be 
impractical to observe this behavior numerically. 



FIG. 5. (Color online) Triplet-singlet mass ratio r as a func- 
tion of 8 and a. (a) The case of S = 1/2. The circle, triangle, 
and down-pointing triangle represent numerically obtained r 
for S = 0.005, 0.01 and 0.015, respectively. The solid, dashed 
and dashed-dotted lines are Eq. (^) with C\ = 0.3. Ms is a 
function of <5 and a, and Ms for a = is used here, (b) The 
case of S — 1. The circle, triangle and down-pointing triangle 
represent numerically obtained r for 8 — S c = —0.005, 0.005 
and 0.01, respectively. S c is determined from Fig. ^j. The solid 
line is Eq. (g) with M s = 0.1J and Ci = 0.6. 



eq. (||), the mass ratio r is 

12V3 

_ V ° + ln(J/M s ) + l/(Ci(Q c -a)) 



1 



±~/3 

ln(.7/Af s ) + l/(Ci(a c -a)) 



(7) 



A fitting of the numerical results with the function (Q) 
is shown in Fig. [s]. The only fitting parameter is the 
non-universal constant C\. For 5=1/2 chain, we use an 
excitation gap with a = as the value of Mg since the 
value of Mg can be estimated through M — u/£, where 
u = nJa/2. The solid, dashed and dashed-dotted lines in 
Fig |(a) are Eq. (@) with d = 0.3 for S = 0.005, 0.01 and 
0.015, respectively. The variation of Eq. (^) by changing 
6 is quite small since the only 5-dependent variable is Ms 
and it is present only inside a logarithm. It is difficult to 
estimate Mg with good precision for 5 = 1 because the 
value of u is not known. However, as we have discussed, 
the Mg dependence is rather weak in Eq. (^|). Thus, in 
a practical range to compare with the numerical results, 
we can set Mg/J = 0.1. Equation ^ with C\ — 0.6 is 
shown as a solid line in Fig |^ (b) . The fitting curves agree 
well with numerical data for both 5 = 1/2 and 5 = 1, in 
the vicinity of a — a c , where the marginal perturbation 
is small. The deviation away from the theory (|7|) can be 
attributed to higher-order correction in both FFPT and 
the renormalization equation. 

Let us come back to the argument in Ref. [17]. As 
we have seen, their idea that r evolves from y3 as 6 
is changed from n mod 27r, does not seem to agree with 
the numerical results. On the other hand, however, where 
the dimerization transition is second order (a < a c ), the 
marginal operator is marginally irrelevant. Thus, in the 
limit 9 — > 7r mod 2-7T (5 — > S c in our spin-chain model), 
the SG theory without the marginal operator becomes 
exact, and r = y/3 should follow. In this sense, their idea 
is still qualitatively correct. However, the marginally 



VII. CONCLUSION 

We have investigated the excitation spectrum of 5 = 
1/2 and 1 frustrated HAF chains with dimerization 5. 
To evaluate particle mass M = u/£, we calculate the 
corresponding correlation function numerically and ex- 
tract the correlation length by using a fitting function 
Ce~ 1 ^ j\fl for a range of large enough and even The 
ratio r of the singlet (the second breather) to the triplet 
(soliton, antisoliton and the first breather) is expected to 
be -\/3 from bosonized SG effective field theory, but r is 
subject to correction from a marginal term, r = \/3~ 
is recovered at the critical next-nearest-neighbor cou- 
pling a = a c , for which the marginal term vanishes. 
At a = a c , the dimerization transition with varying S 
changes from second order, with the critical behavior de- 
scribed by c = 1 CFT, to first order. We give 5 and 
a dependences of r in Eq. (Q) through FFPT and RG 
analysis, r obtained by the iTEBD method is well fitted 
by Eq. (]?]). Our analysis indicates that, for a < a c , the 
mass ratio r asymptotically approaches v3 when i5 — > S c , 
consistently with the argument in Ref. [l7]. However, this 
asymptotic behavior occurs only for exponentially small 
1 5 — <5 C | , and could not be observed in numerical studies 
in the literature and in the present work. 

Finally, we comment on the general-5 case. When the 
dimerization 5 is changed from 1 to — 1 , the phase transi- 
tion happens once for the 5 = 1/2 case (from one dimer- 
ized to the other dimerized phase), and twice for the 
5 = 1 case (from one dimerized to the Haldane phase 
and from the Haldane to the other dimerized phase). In 
the general-5 case, there are 25 transitions from one fully 
dimerized to the other fully dimerized phase, and they are 
the transitions between the partially dimerized phases. 
Around those transition points, the system is represented 
by the same effective field theory as explained in this pa- 
per. 
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